drop if 一套成本==0 & 二套成本==0 & 三套成本==0 & 四套成本==0 & 五套成本==0 & 六套成本==0
drop if 一套面积==. & 二套面积==. & 三套面积==. & 四套面积==. & 五套面积==. & 六套面积==.
drop if 六套面积==.n

replace 六套现值=0 if 六套现值==.
replace 六套税率=0 if 六套税率==.
replace 六套现值=0 if 六套现值==.
replace 减免面积=0 if 减免面积==.
replace 其他市值=0 if 其他市值==.

replace 总税额=一套税率*一套现值+二套税率*二套现值+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率

drop if 一套面积==0 & 二套面积==0 & 三套面积==0 & 四套面积==0 & 五套面积==0 & 六套面积!=0
drop if 二套面积==0 & 三套面积==0 & 四套面积==0 & 五套面积==0 & 六套面积!=0
drop if 三套面积==0 & 四套面积==0 & 五套面积==0 & 六套面积!=0

drop if 一套现值==0 & 二套现值==0  & 三套现值==0 & 四套现值==0  & 五套现值==0 & 六套现值!=0
drop if 二套现值==0  & 三套现值==0 & 四套现值==0 & 五套现值==0 & 六套现值!=0 
drop if  三套现值==0 & 四套现值==0 & 五套现值==0 & 六套现值!=0 
drop if 四套现值==0 & 五套现值==0 & 六套现值!=0
drop if 五套现值==0 & 六套现值!=0

************首套免征***************
replace 首套免征=二套税率*二套现值+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
sum 首套免征

gen 首套纳税能力=(总收入-衣-食-行-总成本-总房贷-首套免征)/首套免征
sum 首套纳税能力
************人均面积减免40***************
replace 人均面积=40*number-一套面积
replace 一套=. if 一套>0  //一套
replace 一套=((-一套)/一套面积)*一套现值*一套税率+二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均面积=. if 人均面积<0 | 人均面积==0
replace 人均面积=人均面积-二套面积
replace 二套=. if 二套>0  //二套
replace 二套=((-二套)/二套面积)*二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均面积=. if 人均面积<0 | 人均面积==0
replace 人均面积=人均面积-三套面积
replace 三套=. if 三套>0  //三套
replace 三套=((-三套)/三套面积)*三套现值*三套税率+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均面积=. if 人均面积<0 | 人均面积==0
replace 人均面积=人均面积-四套面积
replace 四套=. if 四套>0  //四套
replace 四套=((-四套)/四套面积)*四套现值*四套税率+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均面积=. if 人均面积<0 | 人均面积==0

replace 人均面积=0 if 人均面积>0 & 人均面积!=.
replace 三套=四套 if 四套!=.
replace 二套=三套 if 三套!=.
replace 一套=二套 if 二套!=.
replace 人均面积=一套 if 一套!=.
sum 人均面积

************首套免征加人均面积减免40***************
replace 首套面积=40*number-一套面积
replace 一套=. if 一套>0
replace 一套=二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率 if 一套<0 | 一套==0
replace 首套面积=. if 首套面积<0 | 首套面积==0
replace 首套面积=首套面积-二套面积
replace 二套=. if 二套>0
replace 二套=((-二套)/二套面积)*二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率 
replace 首套面积=. if 首套面积<0 | 首套面积==0
replace 首套面积=首套面积-三套面积
replace 三套=. if 三套>0
replace 三套=((-三套)/三套面积)*三套现值*三套税率+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 首套面积=. if 首套面积<0 | 首套面积==0
replace 首套面积=首套面积-四套面积
replace 四套=. if 四套>0
replace 四套=((-四套)/四套面积)*四套现值*四套税率+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 首套面积=. if 首套面积<0 | 首套面积==0

replace 首套面积=0 if 首套面积>0 & 首套面积!=.
replace 三套=四套 if 四套!=.
replace 二套=三套 if 三套!=.
replace 一套=二套 if 二套!=.
replace 首套面积=一套 if 一套!=.

sum 首套面积

************人均价值减免***************
replace 房价=33893 if 省份1==11
replace 房价=15137 if 省份1==12
replace 房价=7039  if 省份1==13
replace 房价=5457  if 省份1==14
replace 房价=6458  if 省份1==21
replace 房价=5748  if 省份1==22
replace 房价=6073  if 省份1==23
replace 房价=24859  if 省份1==31
replace 房价=9070  if 省份1==32
replace 房价=13430  if 省份1==33
replace 房价=6137  if 省份1==34
replace 房价=9284  if 省份1==35
replace 房价=5800  if 省份1==36
replace 房价=6153  if 省份1==37
replace 房价=5038  if 省份1==41
replace 房价=7307  if 省份1==42
replace 房价=4846  if 省份1==43
replace 房价=11416  if 省份1==44
replace 房价=5623  if 省份1==45
replace 房价=11381  if 省份1==46
replace 房价=6605  if 省份1==50
replace 房价=5888  if 省份1==51
replace 房价=4165  if 省份1==52
replace 房价=5664  if 省份1==53
replace 房价=5326   if 省份1==62
replace 房价=5298   if 省份1==63
replace 房价=4243   if 省份1==64

replace 人均价值=房价*number*40-一套现值
replace 一套=. if 一套>0
replace 一套=(-一套)*一套税率+二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均价值=. if 人均价值<0 | 人均价值==0
replace 人均价值=人均价值-二套现值
replace 二套=. if 二套>0
replace 二套=(-二套)*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均价值=. if 人均价值<0  | 人均价值==0
replace 人均价值=人均价值-三套现值
replace 三套=. if 三套>0 
replace 三套=(-三套)*三套税率+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均价值=. if 人均价值<0  | 人均价值==0
replace 人均价值=人均价值-四套现值
replace 四套=. if 四套>0
replace 四套=(-四套)*四套税率+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 人均价值=. if 人均价值<0  | 人均价值==0

replace 人均价值=0 if 人均价值>0 & 人均价值!=.
replace 三套=四套 if 四套!=.
replace 二套=三套 if 三套!=.
replace 一套=二套 if 二套!=.
replace 人均价值=一套 if 一套!=.

sum 人均价值
************人均价值加房贷***************
replace 房贷价值=房价*number*40+总房贷-一套现值
replace 一套=. if 一套>0
replace 一套=(-一套)*一套税率+二套现值*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 房贷价值=. if 房贷价值<0 | 房贷价值==0
replace 房贷价值=房贷价值-二套现值
replace 二套=. if 二套>0
replace 二套=(-二套)*二套税率+三套税率*三套现值+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 房贷价值=. if 房贷价值<0  | 房贷价值==0
replace 房贷价值=房贷价值-三套现值
replace 三套=. if 三套>0 
replace 三套=(-三套)*三套税率+四套税率*四套现值+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 房贷价值=. if 房贷价值<0  | 房贷价值==0
replace 房贷价值=房贷价值-四套现值
replace 四套=. if 四套>0
replace 四套=(-四套)*四套税率+五套税率*五套现值+六套税率*六套现值+其他市值*税率
replace 房贷价值=. if 房贷价值<0  | 房贷价值==0


replace 房贷价值=0 if 房贷价值>0 & 房贷价值!=.
replace 三套=四套 if 四套!=.
replace 二套=三套 if 三套!=.
replace 一套=二套 if 二套!=.
replace 房贷价值=一套 if 一套!=.

sum 房贷价值


*************税额比较***********
egen sum_总税额=sum(总税额)   //总和 
egen sum_首套免征=sum(首套免征) 
egen sum_人均面积=sum(人均面积) 
egen sum_首套面积=sum(首套面积) 
egen sum_人均价值=sum(人均价值) 
egen sum_房贷价值=sum(房贷价值) 


replace 份额1= sum_房贷价值/sum_总税额

*************税负分布***********
egen sum_总收入=sum(总收入)
gen 收入份额=总收入/sum_总收入
sort 收入份额
gen sum_收入份额 = sum(收入份额)   //累加求和
gen 收入阶级=1
replace 收入阶级=2 if sum_收入份额>0.25 & sum_收入份额<=0.5
replace 收入阶级=3 if sum_收入份额>0.5 & sum_收入份额<=0.75
replace 收入阶级=4 if sum_收入份额>0.75 

gen sum_累计收入=sum(总收入)
gen 累计百分比=sum_累计收入/sum_总收入
 
tab 收入阶级 
egen sum_房贷价值=sum(房贷价值)

egen sum_房贷价值1=sum(房贷价值) if 收入阶级==1
gen 税负1=sum_房贷价值1/sum_房贷价值 if 收入阶级==1

egen sum_房贷价值2=sum(房贷价值) if 收入阶级==2
gen 税负2=sum_房贷价值2/sum_房贷价值 if 收入阶级==2

egen sum_房贷价值3=sum(房贷价值) if 收入阶级==3
gen 税负3=sum_房贷价值3/sum_房贷价值 if 收入阶级==3

egen sum_房贷价值4=sum(房贷价值) if 收入阶级==4
gen 税负4=sum_房贷价值4/sum_房贷价值 if 收入阶级==4

sum 税负 if 收入阶级==3
drop if 税负>0.0008 & 收入阶级==2

*************洛伦兹曲线***********
egen sum_总税额=sum(总税额)
gen sum_累计总税额=sum(总税额)
gen 累计税负1=sum_累计总税额/sum_总税额

egen sum_首套免征=sum(首套免征)
gen sum_累计首套免征=sum(首套免征)
gen 累计税负2=sum_累计首套免征/sum_首套免征

egen sum_人均面积=sum(人均面积)
gen sum_累计人均面积=sum(人均面积)
gen 累计税负3=sum_累计人均面积/sum_人均面积

egen sum_首套面积=sum(首套面积)
gen sum_累计首套面积=sum(首套面积)
gen 累计税负4=sum_累计首套面积/sum_首套面积

egen sum_人均价值=sum(人均价值)
gen sum_累计人均价值=sum(人均价值)
gen 累计税负5=sum_累计人均价值/sum_人均价值

egen sum_房贷价值=sum(房贷价值)
gen sum_累计房贷价值=sum(房贷价值)
gen 累计税负6=sum_累计房贷价值/sum_房贷价值

twoway line 累计税负1 sum_收入份额 , sort || line 累计税负2 sum_收入份额 , || line 累计税负3 sum_收入份额 ,|| line 累计税负4 sum_收入份额 , || line 累计税负5 sum_收入份额 , || line 累计税负6 sum_收入份额 , || line sum_收入份额 sum_收入份额 , ///
xlabel(0(0.1)1) ylabel(0(0.1)1)      ///
xline(0(0.1)1) yline(0(0.1)1)        ///
legend(label(1 "洛伦兹曲线") label(2 "绝对平等线")) ///
plotregion(margin(zero)) aspectratio(1) scheme(s1mono)

sum 税负 if 收入阶级==1
drop if 税负==0 & 收入阶级==3
drop if 税负>0.00015 & 收入阶级==1





